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We discuss the influence of fluctuations on thermodynamics near the chiral phase tran- 
sition within Polyakov loop extended quark-meson model based on the functional renor- 
. malization group (FRG) method. We include the gluon fields in the FRG flow equation 

self-consistently on the mean-field level. We focus on the properties of the phase diagram 
and net-baryon number fluctuations. 
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§1. Introduction 



Understanding the phase structure and critical properties of strongly interacting 
matter is one of the central problems addressed in studies of QCD thermodynamics. 
In the context of heavy ion experiments of particular importance is to find measur- 
able observables that can be used as probes of the deconfinement and chiral phase 
transitions expected in QCDP Effective chiral models were shown to be very useful 
to study critical phenomena in the strongly interacting system and to quantify the 
phase transition due to deconfinement and chiral symmetry breaking!^ 

Particularly interesting are models like the Polyakov loop extended Nambu- 
Jona-Lasinio (PNJL)PJ or quark-meson (PQM) 4 -' models which can account for both 
£S) \ deconfinement and chiral symmetry breaking. Both models can reproduce essential 

properties of QCD thermodynamics obtained in the Lattice Gauge Theory already 
within the mean-field (MF) approximation. However, to correctly account for criti- 
cal behavior and scaling properties near the chiral phase transitions one needs to go 
beyond the mean-field approximation and include fluctuations and non-perturbative 
dynamics. This can be achieved by using methods based on the functional renor- 
malization group (FRG) I 

In this paper we summarized the influence of non-perturbative effects in the 
presence of a gluonic background on the thermodynamics the near chiral phase tran- 
sition based on the PQM model. We compare the phase diagram and properties 
of net-quark number fluctuations obtained within the FRG and MF calculations. 
Such a comparison quantifies explicitly the influence of mesonic fluctuations near 
the chiral phase transition at finite temperature and density.^'® 

§2. The FRG flow equation in a background gluon field 

The role of non-perturbative effects in the thermodynamic description of the 
PQM model can be studied within the functional renormalization group (FRG) 
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Fig. 1. The phase diagrams for the PQM model in the mean-field approximation (left panel) 
and in the functional renormalization group approach (right panel). The arrows indicate lines 
corresponding to different values of \xjT. The shaded lines define a critical region (see text). 



method. The FRG is based on an infrared regularization with the momentum scale 
parameter of the full propagator which turns the corresponding effective action into 
the scale fc-dependent functional r k ® In general, in the PQM model, the formula- 
tion of the FRG flow equation requires an implementation of the Polyakov loop as a 
dynamical field. However, in the current formulation we treat the Polyakov loop as a 
background field which is introduced self-consistently on the mean-field level!^ Fol- 
lowing our previous worl^® we formulate the flow equation for the scale-dependent 
grand canonical potential as 



d k f2 k (£,t;T,fj,) 




— l + 2n B (E 7r ;T) 



1 



l + 2n B (E a ;T) 



(2-1) 



l-N(e,f,T,n)-N{£,t;T,n) 



Here n B (E 7Tta ; T) is the bosonic distribution function 

1 



n B (E ni(7 ;T) 



exp(E Vta /T) - 1 



with the pion, E n = y k 2 + f2 k , and sigma, E a = y k 2 + Q k + 2p Q k energies, where 

the primes denote derivatives with respect to p = (a 2 + tt 2 )/2 and 17 = fl + ca with 
c = m 2 /^ being the external symmetry breaking term. The functions N (£,£*; T, p) 
and N(£, £*; T, p) which are defined by 



N(£,t;T,p) 
N(£,t;T,p,) 



1 + 21* exp[P(E q - //)] + eexp[2P(E q - p)] 



l + 3£exp[2(3(E q 
N(t,£;T,-p), 



n)] + Wexp\p(E q 



(2-2) 



are fermionic distributions which are modified because of coupling to gluons. The 
E q = \Jk 2 + 2g 2 p is the quark energy. The minimum of the thermodynamic potential 
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is determined by the stationarity condition 



da 



dflk 



da 



- c = 0. (2-3) 

<7=C7 fc 



The flow equation (|2-ip is solved numerically with the initial cutoff A = 1.2 GeV 
( see details in RefP®). The initial conditions for the flow are chosen to reproduce 
the physical pion mass m n = 138 MeV, the pion decay constant f n = 93 MeV, the 
sigma mass m a = 600 MeV and the constituent quark mass m q = 300 MeV at the 
scale k = and for zero temperature and chemical potential T = fj, = 0. 

By solving equation (|2-ip one obtains the thermodynamic potential for quarks 
and mesons, f2k-to(£, £*'■> T, (J,), as a function of Polyakov loop variables, £ and £*. The 
full thermodynamic potential f2(£, £*;T, //) in the PQM model which includes quarks, 
mesons and gluons degrees of freedom is obtained by adding to f2k->o(£, £*;T, fj,) the 
effective gluon potential, U(£,£*)/T 4 = -b 2 (T)£*£ - b 3 (£ 3 + £* 3 ) + b 4 (£*£) 2 . Thus, 

{}(£,£* -T^) = Q k ^ Q {t,t-T,n) +U(£,£*), (2-4) 

with parameters b{ chosen to reproduce the equation of state of the pure SU C (3) 
lattice gauge theoryP 

At a given temperature and chemical potential, the Polyakov loop variables, £ 
and £*, are determined by the stationarity conditions: 

^- £ f2(£,£*;T,n) = 0, ^f2(£,t;T, n) = 0. (2-5) 



The thermodynamic potential (|2-4p does not contain contributions of statistical 
modes with momenta larger than the cutoff A. In order to obtain the correct high- 
temperature behavior of thermodynamics we need to supplement the FRG potential 
with the contribution of high-momentum states. For k > A we use the flow equation 
for quarks interacting with the Polyakov loops P 

d k f2£(T, M ) = ^fjf [l - N(£, t;T, n) - N(£, f ; T, /i)] (2-6) 

where the dynamical quark mass is neglected. 

To obtain the complete thermodynamic potential of the PQM model we integrate 
Eq. (|2-6p from k = oo to k = A where we switch to the PQM flow equation (|2-ip . 



§3. Thermodynamics of the PQM model in the presence of quantum 
fluctuations and gluonic background 

In order to illustrate the importance of mesonic fluctuations and gluonic back- 
ground on thermodynamics in the PQM model we will compare the FRG results 
with those obtained under a mean-field approximation. In the latter, both quantum 
and thermal fluctuations are neglected and the mesonic fields are replaced by their 
classical expectation values. The formulation of the thermodynamic potential in the 
MF approximation can be found in RefsP® 
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Fig. 2. The coefficient C2 as a function of temperature for different values of fi/T for the PQM 
model in the mean-field approximation (left panel) and in the FRG approach (right panel). 
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Fig. 3. The coefficient C4 as a function of temperature for different values of fi/T for the PQM 
model in the mean-field approximation (left panel) and in the FRG approach (right panel). 



Fig. Q] shows the phase diagrams of the PQM model obtained in the FRG ap- 
proach and in the mean-field approximation. For the physical pion mass and mod- 
erate values of the chemical potential the PGM model exhibits a smooth cross-over 
chiral transition. In Fig. [I] we define the transition region as a band where the tem- 
perature derivative of the order parameter exhibits 5%-deviations from its maximal 
value. At larger \x the cross-over line terminates at the CEP where the transition is 
second-order and belongs to the universality class of the three dimensional Z(2)-spin 
system. 

Comparing the resulting phase diagrams of the PQM model obtained within MF 
and FRG approach we find a clear shift of the position of the chiral phase boundary 
to higher temperatures due to mesonic fluctuations. Changes in the phase diagram 
were previously reported in the quark-meson model within the FRG approach.® 
However, here due to gluonic background, which is explicitly included in our FRG 
calculations, we find a significant shift of the CEP to higher temperature.® 

The fluctuations of conserved charges were shown to be an excellent probe of the 
chiral and deconfinement transition expected in strongly interacting mediumP'EB 
A particular role was attributed to net-baryon number fluctuations and their higher 
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Fig. 4. The kurtosis i?4,2 as a function of temperature for different /j,/T calculated in the PQM 
model under the mean-field approximation (left panel) and in the FRG approach (right panel). 



moments which are characterized by the generalized susceptibilities, 

c n(T) = Wr)B . (3-1) 

The temperature and chemical potential dependence of the net-quark number suscep- 
tibility C2 in the PQM model is shown in Fig.[2]for the MF and the FRG calculations. 
At vanishing \i the coefficient C2 increases monotonously with temperature, similar to 
the thermodynamic pressure. However, at finite chemical potential, the susceptibil- 
ity C2 develops a peak structure. The amplitude of this peak increases with chemical 
potential towards the CEP where C2 diverges. In the high temperature phase the C2 
converges to the Stefan Boltzmann value, cf B = 2[1 + 2(3/i/-7rT) 2 ]. 

As it is seen in Fig. [2] the peak structure in C2 is more pronounced in the mean- 
field at n/T = 1 than in the FRG at fi/T =1.5. This is in spite of the fact that the 
location of the CEP in the FRG is closer to the fj,/T = 1.5 than the corresponding 
one for the mean-field approximation at fi/T = 1 (see Fig. [1]). This shows that 
the criticality of C2 as a function of a distance to the CEP appears earlier in the 
mean-field than in the FRG calculations. 

The generalized susceptibility C2 is sensitive to changes in the chemical potential 
and is clearly influenced by meson fluctuations and the gluon background. The C2 is 
strictly positive for all values of // and T. At finite [i the positivity is not preserved 
for c n -moments with n > 2. 

Fig. E] shows the fourth order cumulant moment of the net-quark number for 
different values of fJ,/T. The fourth order cumulant is strictly positive for vanish- 
ing chemical potential. However, for higher values of (jl/T, C4 becomes negative 
in the vicinity to the cross-over transition. The chemical potential independent 
Stefan-Boltzmann limit cf B = 2N C N f /it 2 is reproduced at temperatures T » T c . 
Comparing the MF with the FRG results for C4 it is clear that meson fluctuations 
are essentially modifying properties of quark susceptibilities. In the transition re- 
gion C4 is suppressed in the FRG relative to the MF results. The meson fluctuations 
quantified by the FRG method provide smoothing of net-quark susceptibilities near 
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the crossover transition which is substantially broader in the FRG than in the MF 
approach. 

Discussing fluctuations of conserved charges, a particular role is attributed to 
the so called kurtosis of net-quark number fluctuations which is given by the ratio, 
i?4 ; 2 = 04/02-^ This key observable is not only sensitive to the chiral but also to 
the confinement-deconfinement transition. In the asymptotic regime of high and 
low temperature and at vanishing chemical potential the kurtosis reflects the quark 
content of the baryon-number carrying effective degrees of freedom.^ Therefore, at 
low T in the confined phase, R^2 — = 9 while in an ideal gas of quarks R^2 — 1- 
This property was shown to be independent from the particular form of the mass 
spectrum and the number of degrees of freedom in a system. 

Fig. H] shows the C4/C2 ratio calculated as a function of temperature along 
different paths in the (T,fj,)— plane quantified by fixed n/T. In the PQM model the 
kurtosis shows an expected drop from — 9 to R^ — 1 in the transition region. 
At vanishing \i it exhibits a peak at the transition temperature. The height of this 
peak depends not only on the pion mas^ 1 ® but is also influenced by the value of 
the chemical potential. The mesonic fluctuations weaken the peak structure both 
at finite and at vanishing quark density. For finite \i the kurtosis becomes negative 
following the same trends as seen in the fourth order cumulant. 

Summary 

We have discussed the influence of quantum fluctuations near the chiral phase 
transition within the Polyakov loop extended quark-meson (PQM) model based 
on the functional renormalization group method (FRG). We have shown that non- 
perturbative dynamics introduced in the FRG approach essentially modifies predic- 
tions of the PQM model derived in the mean-field approximation. In particular, we 
have demonstrated quantitative changes of the phase diagram and modification of 
net-quark number fluctuations and their higher moments. We have indicated the 
role and importance of the ratio of the fourth to second order cumulant moments to 
identify deconfinement and chiral phase transition in a strongly interacting medium. 

References 

1) M. M. Agga rwal et al. [STA R Collaboration], I arXiv: 1007.26131 [nucl-ex]. F. Karsch and 
K. Redlich. larXiv7l 007.2581 [hep-ph]. 

2) C. Sasaki, Nucl. Phys. A 830 (2009) 649C and references therein. 

3) K. Fukushima, Phys. Lett. B 591, 277 (2004). C. Ratti, M. A. Thaler and W. Weise, Phys. 
Rev. D 73, 014019 (2006). C. Sasaki, B. Friman and K. Redlich, Phys. Rev. D 77, 034024 
(2008). Phys. Rev. Lett. 99, 232301 (2007). Phys. Rev. D 75, 074013 (2007). 

4) B. J. Schaefer, J. M. Pawlowski and J. Wambach, Phys. Rev. D 76, 074023 (2007). 
V. Skokov, et al., Phys. Rev. D 82 (2010) 034029. 

5) For a review see e.g., J. Berges, N. Tetradis and C. Wetterich, Phys. Rept. 363, 223 (2002) 

6) B. Stokic, B. Friman and K. Redlich, Eur. Phys. J. C 67 (2010) 425. E. Nakano, B. J. Schae- 
fer, B. Stokic, B. Friman and K. Redlich, Phys. Lett. B 682 (2010) 401. 

7) V. Skokov, B. Stokic, B. Friman and K. Redli ch, Phys. Rev. C 82 (2010) 015206 

8) V. Skokov, B. Friman and K. Redlich, larXiv: 1008.4 570 [hep-ph]. 

9) T. K. Herbst, J. M. Pawlowski and B. J. Schaefer, arXiv: 1008.0081 [hep-ph]. 
10) S. Ejiri, F. Karsch and K. Redlich, Phys. Lett. B 633, 275 (2006). 



